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ABSTRACT 


The  range  of  fluid  dynamics  applications  in  the  lung 
is  vast,  encompassing  many  interesting  problems  rarely 
touched  in  any  other  single  engineering  area.   A  survey 
of  recent  work  in  this  field  was  conducted.   Additionally, 
a  numerical  solution  for  gas  mixing  in  a  simple  model  of  the 
lung  was  developed  indicating  its  application  to  a  wide 
range  of  similar  problems.   Finally,  new  parameters  were 
experimentally  investigated  in  flow  through  a  thin-walled 
collapsible  tube.   Results  indicated  the  effects  of  the 
parameters  to  be  significant  in  the  flow  of  blood  in  the 
pulmonary  capillaries  as  controlled  by  the  surrounding 
aveolar  pressure. 
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I.   SURVEY  OF  FLUID  DYNAMICS  IN  THE  LUNG 

A.  INTRODUCTION 

The  study  of  fluid  dynamics  in  the  lung  is  a  wide  field, 
encompassing  a  large  range  of  problems.   The  vastness  of  the 
field  is  easily  understood  when  the  complexity  of  the  organi- 
zation of  the  respiratory  tract  and  the  cardiovascular  system 
is  considered.   Many  of  the  problems  which  arise  in  these 
areas  are  unique  even  when  the  fluid  dynamics  of  the  problems 
are  considered  in  isolation.   Often,  these  problems  require 
work  in  areas  of  fluid  dynamics  that  have  never  been 
attempted  in  purely  engineering  applications.   Frequently, 
new  researches  in  the  field  center  on  very  basic  fluid  flow 
problems,  and  the  results  are  as  interesting  as  any  that  are 
suggested  by  engineering  applications.   Figure  1  indicates 
the  main  areas  of  study  within  the  cardiovascular-respiratory 
system. 

B.  FLUID  FLOW  IN  LARGE  AIRWAYS  AND  ARTERIES 

In  both  the  large  airways  of  the  lung  and  in  the  larger 
arteries  the  Reynolds  number  may  vary  from  several  thousand 
to  one  or  two  hundred.   This  range  of  flows  has  not  been 
closely  studied,  particularly  the  large  Reynolds  number  end. 
In  most  engineering  applications  large  Reynolds  number 
studies  have  been  limited  to  fully  developed  turbulent  flow. 
In  the  lung,  however,  high  Reynolds  number  flows  are  only 
sporadically  or  locally  turbulent,  seldom  becoming  fully 
developed. 


In  studies  by  Kuchar  and  Ostrach  [Ref.  1]  for  large 
Reynolds  number  flow  regimes  including  the  effect  of 
viscosity,  pulsatile  flow,  and  vessel  elasticity,  the 
developing  flow  was  found  to  have  a  strong  boundary  layer 
character,  resulting  in  high  wall  shear  stresses  near  the 
vessel  entrances.   The  development  length  was  found  to  be 
large  and  dependent  on  the  Reynolds  number  and  unsteadiness 
parameter,  while  wall  distensibility  was  not  significant. 
Thus  flow  development  characteristics  are  significant  in  the 
large  airways  and  arteries. 

Additionally,  the  large  vessels  may  be  subject  to 
curvature,  inducing  large  secondary  flows  due  to  centrifugal 
forces.   These  forces,  acting  as  a  result  of  curvature  of 
the  primary  flow  streamlines,  produce  secondary  flow 
velocities  at  right  angles  to  the  primary  streamlines. 
Study  by  McConologue  and  Srivastava  [Ref.  2]  has  indicated 
that  the  resulting  motion  of  a  central  fluid  is  toward  the 
outside  of  the  bend,  and  a  return  flow  toward  the  inside  of 
the  bend  near  the  wall.   The  motion  toward  the  outside  wall 
is  nearly  uniform  while  the  return  flow,  retarded  by  viscous 
stress,  is  limited  to  a  thin  boundary  layer.   A  result  of 
this  fact  was  the  increase  in  critical  Reynolds  number 
based  on  pipe  diameter  for  transition  in  curved  pipes 
compared  with  straight  pipes.   In  addition,  the  secondary 
flow  promotes  a  much  more  rapid  lateral  mixing  than  would 
be  possible  for  normal  laminar  flow. 


Work  has  been  undertaken  to  study  the  combined  effect 
of  the  last  two  topics  due  to  the  fact  that  the  first  entry 
region  to  both  the  cardiovascular  system  and  respiratory 
tract  is  a  curved  tube:   the  great  arch  of  the  Aorta,  and 
the  mouth  and  trachea. 

C.   COMPLEX  BRANCHING 

Complex  branching  of  the  passages  exists  in  both  the 
respiratory  tract  and  the  cardiovascular  system.   As  the 
flow  branches  early  in  each  system  it  undergoes  a  large 
distortion  at  each  branch,  and  the  flow  may  not  recover 
before  reaching  the  next  bifurcation.   In  the  smaller  air- 
ways and  capillaries  where  the  Reynolds  number  is  too  small 
for  this  effect,  the  flow  becomes  difficult  to  analyze 
due  to  the  large  number  of  tubes  in  parallel  and  the  lack 
of  information  on  the  distribution  of  tubes  by  length  and 
diameter. 

Pedley,  Schroter,  and  Sudlow  [Ref.  3]  found  that  in 
the  higher  Reynolds  number  airways  secondary  flow  was 
exhibited  downstream  of  each  bifurcation,  with  peak  axial 
velocities  just  outside  the  boundary  layer  on  the  inner 
wall  of  the  junction.   These  secondary  motions  decayed 
slowly  downstream. 

Additional  study  by  Caro,  Fitz-Gerald,  and  Schroter 
[Ref.  4]  has  indicated  that  the  rate  of  shear  in  a  branching 
system  was  dependent  on  two  factors:   a  flow  factor  depend- 
ing on  the  total  flow  in  a  particular  tube,  and  a  geometrical 


factor.   The  geometrical  factor  was  found  to  decrease 
downstream.   It  was  also  greater  on  the  outside  of  a  bend 
than  on  the  inside.   At  a  bifurcation  the  factor  was  largest 
on  the  central  flow  divider,  the  outer  bends  of  the  two 
daughter  branches  becoming  inside  bends,  effectively.   The 
flow  factor  was  also  found  to  decrease  with  distance  down- 
stream.  The  general  result  was  that  the  rate  of  shear  was 
least  nearer  the  heart  or  mouth  and,  within  these,  at  the 
site  where  the  geometrical  factor  was  smallest. 

In  steady  flow  a  branching  system  affects  pressure 
drops  at  each  series  of  bifurcations.   Here  the  effect  of 
resistance  differs  for  the  arterial  problem  and  the 
bronchial  problem,  as  discussed  by  Lighthill  [Ref .  5] .   In 
the  bloodstream  the  greatest  pressure  drops  occur  in  the 
smaller  vessels  of  diameter  less  than  0.1  mm.,  whereas  in 
the  airways  the  maximum  pressure  drops  occur  early  in  the 
branching,  soon  after  the  fifth  generation.   The  airways 
maintain  a  much  larger  daughter-to-mother  branch  ratio 
closer  to  the  system  entrance.   This  coupled  with  the  higher 
Reynolds  number  in  those  regions  would  tend  to  cause  separa- 
tion, a  possible  explanation  of  the  phenomenon.   The  arterial 
resistance  problem  is  said  to  be  a  low  Reynolds  number 
problem  while  bronchial  resistance  is  considered  a  high 
Reynolds  number  problem. 

D.   BLOOD  MICROCIRCULATION 

Flow  in  the  small  capillaries  of  the  lung  differs 
greatly  from  that  usually  found  in  engineering  fluid 


dynamics.   The  feature  which  distinguishes  this  type  of 
flow  is  that  the  fluid  inertial  effects  are  totally  negli- 
gible.  The  Reynolds  number  is  so  low  that  the  flow  pattern 
is  determined  solely  by  the  need  for  pressure  gradients  and 
viscous  forces  to  be  in  exact  equilibrium.   Under  these 
conditions,  Bernoulli's  equation  is  not  applicable  and  there 
is  essentially  no  difference  between  static  and  dynamic 
pressures.   Centrifugal  forces  are  also  negligible:   fluid 
can  negotiate  sharp  bends  with  no  secondary  flow  effects. 
In  general,  fluid  motion  is  insensitive  to  vessel  geometry 
and  there  is  little  tendency  for  separation.   This  insensi- 
tivity  to  geometry  was  further  demonstrated  by  Fitz-Gerald 
[Ref.  6]  with  respect  to  entrance  length  effects.   The 
difference  between  Poiseuille  motion  of  a  fluid  and  fluid 
entering  a  tube  at  a  large  angle,  with  the  same  total  flow, 
was  reduced  by  an  order  of  magnitude  in  an  axial  distance 
of  half  the  tube  radius. 

Propagation  of  changes  in  pressure,  in  the  large  vessels 
governed  mainly  by  a  balance  between  distensibility  and 
inertia,  is  governed  in  the  smaller  vessels  entirely  by  a 
balance  between  distensibility  and  resistivity.   A  feature 
of  the  blood  circulation  in  most  parts  of  the  body  is  the 
vasomotor  control  of  peripheral  perfusion,  the  ability  of 
smooth  muscle  to  alter  the  diameter  of  small  vessels  in 
order  to  vary  the  flow  rate  through  the  vessels  by  an  order 
of  magnitude.   Contrary  to  the  systemic  microcirculation 
system,  microcirculation  in  the  lung  exhibits  almost  no 
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vasomotor  control   in  the  pre-capillary  vessels  or  elsewhere. 
As  noted  by  Lighthill  [Ref .  5] ,  most  of  the  resistance  to 
flow  appears  to  be  in  the  capillaries  themselves. 

Perfusion  in  the  capillaries  can  vary  greatly  with 
location  in  the  lung,  an  effect  based  predominantly  on 
gravity  as  presented  by  Glaister  [Ref.  '7] .   A  typical 
pulmonary  capillary  is  easily  distorted  by  pressure  dif- 
ferences between  the  capillary  and  alveolus.   In  effect, 
the  capillary  is  a  thin-walled  collapsible  tube  in  which 
flow  is  governed  by  the  alveolar  pressure  P  ,  and  the 
upstream  and  downstream  capillary  pressures,  P.  and  P„, 
respectively.   The  following  information  was  obtained  and 
verified  by  Conrad  [Ref.  8],  Holt  [Ref.  9],  and  many 
others.   The  effect  of  gravity  was  to  divide  the  lung  into 
three  basic  zones.   In  zone  1,  at  the  top  of  the  lung,  P 
was  greater  than  P,  or  P„,  and  there  was  no  flow.   In  zone  2, 
the  middle  of  the  lung,  P~  was  less  than  P  which  was  less 

L,  cl 

than  P, ,  and  flow  was  governed  by  the  pressure  difference 

p, -p  .   In  zone  3,  low  in  the  lung,  P  was  less  than  P,  or 
*±  *a  '  6*   a  1 

P~,  and  flow  was  unobstructed  and  proportional  to  the 
pressure  difference  P-.-P-.   An  experimental  investigation  of 
flow  through  a  thin-walled  collapsible  tube  is  included  in 
this  paper  in  Section  III. 

E.   FLUID  PROPERTIES 

A  problem  relatively  unique  to  physiological  fluid 
flow  lies  in  the  unusual  fluid  properties  of  blood,  and  to 
a  lesser  extent  of  air.   Blood  is  a  suspension  (40-501  by 
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volume)  of  small  deformable  bodies,  mainly  red  cells,  in  a 
transparent  fluid,  the  blood  plasma.   Although  the  flow 
properties  of  plasma  are  similar  to  those  described  by 
Newton's  viscosity  law,  the  concentration  of  red  cells  at 
low  strain  rates  yields  substantially  increasingly  effective 
viscosities.   However,  Lighthill  [Ref.  5]  indicates  that 
data  suggest  that  at  the  higher  strain  rates  in  the  larger 
vessels  the  errors  are  not  large  in  treating  the  blood  as 
a  Newtonian  fluid  with  constant  viscosity.   Studies  of 
particulate  suspensions  by  Mason  and  Goldsmith  [Ref.  10] , 
simulating  blood  over  a  range  of  Reynolds  numbers  extending 
from  creeping  viscous  flow  to  that  where  inertial  effects 
play  an  important  role,  indicated  that  at  low  Reynolds 
numbers  for  a  single  particle,  particle  deformation  and 
interaction  with  the  wall  forced  migration  to  the  axis  of 
the  tube.   At  higher  Reynolds  numbers,  particles  also 
migrated  laterally,  but  to  an  eccentric  radial  position. 
In  blood-similar  suspensions  at  low  Reynolds  numbers, 
blunting  of  the  velocity  profiles  was  evident,  the  dispersions 
were  non-Newtonian,  and  cells  were  deformed  significantly. 
At  higher  Reynolds  numbers  in  steady  flow,  and  in  oscillatory 
and  pulsatile  flow,  cells  tended  to  migrate  inward  from  the 
tube  wall.   This  resulted  in  a  two-phase  flow  of  a  central, 
more  viscous  core  surrounded  by  a  thin,  particle-depleted 
and  less  viscous  zone,  with  an  attendant  decrease  in  the 
energy  dissipated  in  the  tube. 
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Applications  of  "hydrodynamic  lubrication  theory"  by 
Lighthill  [Ref.  11]  to  blood  flow  in  the  narrowest  capil- 
laries have  corresponded  to  the  observed  facts  that  the 
stress  necessary  to  deform  a  capillary  wall  was  nearly  two 
orders  of  magnitude  larger  than  that  required  to  deform  a 
red  cell.   Thus  it  is  that  the  cells  deform  very  substan- 
tially in  order  to  pass  through  the  smallest  capillaries. 
Lubrication  theory  further  suggests  that  the  stress,  which 
the  wall  of  a  very  narrow  capillary  passage  must  transmit 
to  a  red  cell  sliding  past  it  so  that  the  cell  will  deform 
sufficiently,  must  be  transmitted  by  viscous  action  in  a 
thin  lubricating  layer  of  plasma^  that  the  plasma  film 
thickness  must  adjust  passively  to  the  value  required;  and 
that  sliding  may  cease  when  the  plasma  film  thickness 
becomes  too  small.   Theoretical  observations  also  indicated 
that  the  film  thickness  should  vary  with  the  square  root  of 
the  cell  velocity.   This  is  an  interesting  non-linear  effect. 

Air  is  also  a  .suspension  of  dust  and  other  fine 
particles,  and  significant  study  has  been  devoted  to  the 
behavior  of  airborne  particles  in  the  lung.   The  deviation 
of  particle  motion  from  gas  motion  is  caused  by  Brownian 
motion,  gravity,  and  inertial  forces.   Altshuler  [Ref.  12] 
stated  that  for  particles  of  less  than  0.2  ym  the  dominant 
force  is  Brownian  motion.   For  particles  of  greater  than  1  ym 
the  dominant  effects  are  settling  and  inertial  displacement. 

Owen  [Ref.  13]  observed  that  for  inspired  flow  rates  of 
40  L./min.  and  greater  the  boundary  layer  on  the  tracheal 
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wall  was  largely  turbulent.   In  addition,  if  the  air 
entering  the  trachea  separates  at  the  vocal  chords,  almost 
fully  developed  turbulent  pipe  flow  may  exist  and  decay 
slowly.   This  would  suggest  that  many  particles  are  deposited 
on  the  tracheal  walls;  however,  the  normal  rate  of  evapora- 
tion from  the  tracheal  wall  causes  a  radial  motion  sufficient 
to  impede  the  deposition  of  all  but  the  finest  or  coarsest 
particles . 

G.   GAS  MIXING  IN  THE  RESPIRATORY  TRACT 

At  ordinary  breathing  frequencies  the  main  mechanical 
impedance  to  gas  flow  within  the  lungs  stems  from  the  elastic 
distensibility  of  the  air  spaces  terminating  the  airways 
(alveoli) ,  the  flow  resistance  of  the  air  passages  being 
extremely  small.   Due  to  the  weight  of  the  lungs,  there  is 
a  gradient  of  increasing  distension  from  bottom  to  top. 
Mead  [Ref.  14]  found  that  the  lower,  less  distended  regions 
were  more  compliant  than  the  upper,  more  distended  ones,  and 
that  there  was  a  gradient  of  increasing  gas  flow  from  top  to 
bottom.   As  a  result,  the  upper  regions  of  the  lungs  received 
only  half  the  bulk  flow  of  gas  per  gram  of  lung  tissue 
received  by  the  lower  regions. 

Davidson  and  Fitz-Gerald  [Ref.  15]  discovered  that  for 
various  geometrical  models  for  small  airways  in  the  lung, 
the  bulk  flow  was  significant  for  gas  transport  in  the  final 
generations  of  the  lung. 

In  single-breath  nitrogen  washout  experiments,  as 
discussed  by  Farhi  [Ref.  16]  and  Pedley  [Ref.  17],  in  which 
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the  subject  inspired  a  single  breath  of  pure  oxygen,  then 
expired  it,  the  N?    concentration  was  measured  at  the  mouth 
during  expiration.   The  results  indicate  that  initially 
there  is  no  N-  because  the  expired  gas  is  the  pure  0?  most 
recently  inspired,  which  never  made  contact  with  the  N? 
rich  air  already  present  in  the  lung  before  0   inspiration 
began.   Following  this  region,  the  N?  concentration  rapidly 
increased  to  a  nearly  constant  value  on  the  "alveolar 
plateau."   The  alveolar  plateau  represented  air  from  deep 
in  the  lungs  where  the  N~  concentration  was  initially  that 
of  ordinary  air. 

Gas  mixing  in  a  simple  model  of  the  lung  was  solved 
numerically  in  Section  II  of  this  paper. 

G.   GAS  EXCHANGE  AT  THE  ALVEOLAR  INTERFACE 

The  interchange  of  materials  (0~,  CO- )  by  the  circulatory 
system  depends  on  two  processes:   (1)  circulation  of  blood 
through  the  capillary  beds  and  (2)  transport  of  various 
components  of  the  blood  through  the  capillary  walls.   Net 
transport  across  capillary  walls  is  driven  by  two  basic 
forces:   (1)  differences  in  hydrostatic  or  effective 
osmotic  pressure  of  the  fluids  on  either  side  (ultrafiltra- 
tion transport)  and  (2)  differences  in  concentrations  of 
permeating  solutes  on  either  side  (diffusion  transport) . 

Renkin  .[Ref.  18]  found  that  control  of  ultrafiltration 
transport  was  brought  about  by  changes  in  the  ratio  of  the 
resistance  of  the  latter  half  of  a  capillary  bed  to  the 
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total  resistance  of  the  bed,  and  by  changes  in  the  number 
of  open  capillaries.   Control  of  diffusion  transport  was 
brought  about  by  changes  in  the  number  of  open  capillaries, 
and  by  changes  in  local  distribution  of  capillary  blood 
flow  relative  to  capillary  surface  and  permeability. 

As  blood  flows  along  the  pulmonary ' capillaries  it 
receives  oxygen  from  the  alveoli,  and  delivers  carbon 
dioxide  to  them.   Thus  concentration  gradients  for  the  gases 
develop  along  the  capillaries,  increasing  for  oxygen  and 
decreasing  for  carbon  dioxide.   Study  of  digital  computer 
models  of  pulmonary  gas  exchange  by  Milhorn  and  Pulley 
[Ref.  19]  revealed  that  the  partial  pressure  of  blood  reached 
equilibrium  with  alveolar  oxygen  at  about  two-fifths  of  the 
total  distance  along  the  capillary  pathway.   Also  the  partial 
pressure  of  the  CO-  reached  equilibrium  with  the  blood  about 
halfway  along  the  capillary  pathway.   Additionally,  the 
partial-pressure  gradient  of  oxygen  had  a  concave  shape  due 
to  the  convex  shape  of  its  dissociation  curve.   The  concavity 
did  not  exist  in  either  the  concentration  gradient  of  oxygen 
or  in  each  of-  the  curves  for  carbon  dioxide. 

The  effect  of  pulsatile  flow  of  blood  on  oxygenation  of 
the  blood  in  pulmonary  capillaries  for  an  ideal  lung  was 
investigated  by  Graeser,  Kim,  and  Crandall  [Ref.  20]  for 
several  combinations  of  exercise  and  atmospheric  oxygen 
content.   The  pulsatile  flow  was  found  to  give  substantially 
the  same  end  capillary  alveolar-arterial  tension  gradient 
as  that  found  for  a  constant  velocity  model.   When  breathing 
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air,  the  pulsatile  flow  effect  yielded  a  10-15%  greater 
diffusing  capacity  than  the  constant  velocity  model,  but 
no  effect  was  observed  in  low  oxygen  cases. 
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II.   NUMERICAL  SOLUTION  FOR  GAS  MIXING 
IN  A  SIMPLE  MODEL  OF  THE  LUNG 

A.   INTRODUCTION 

The  diffusion  of  gas  throughout  the  lung  is  a  complex 
problem  dependent  on  both  convective  transfer  and  molecular 
diffusion.   A  typical  physiological  experiment  in  this  field 
is  the  single  breath  nitrogen  washout  experiment  described 
by  Farhi  [Ref.  16]  and  Pedley  [Ref.  17].   The  subject, 
starting  from  an  arbitrary  lung  volume,  inspired  a  single 
breath  of  pure  oxygen.   The  subject  then  exhaled  and  the 
N2  concentration  at  the  mouth  was  measured  during  expira- 
tion.  The  resulting  graph  of  N-  concentration  versus  time 
was  typically  similar  to  that  of  Figure  2,  in  which  high- 
frequency  oscillations  due  to  the  cardiac  cycle  have  been 
averaged  out.   In  region  I  there  was  essentially  no  N~  due 
to  the  fact  that  the  expired  gas  was  the  pure  0?  most 
recently  inspired,  which  never  reached  the  N?  rich  air 
already  in  the  lung  upon  inspiration.   Region  II  was  a 
transition  region  where  the  N?  concentration  increased 
rapidly  with  time,  reaching  a  nearly  constant  level  called 
the  "alveolar  plateau"  in  region  III.   The  alveolar 
plateau  represented  air  from  deep  in  the  lungs  where  the 
initial  N?  concentration  was  that  of  ordinary  air.   If  the 
total  volume  of  0?  inspired  (V.  )  (related  to  time  by  the 
mass  flow  rate)  was  less  than  the  anatomical  dead  space  (V~D) 
(the  volume  of  the  airways  as  opposed  to  the  alveoli) ,  then 
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the  volume  V"n  at  which  the  plateau  was  reached  was 
approximately  equal  to  V.  .   If  V.   was  greater  than  V~n, 
Vn  was  approximately  equal  to  Vn,  although  the  exact  relation- 
ship was  a  function  of  the  gas  mixing  and  diffusion  of  the 
inspired  gas  with  the  air  previously  present. 

A  significant  effect  in  region  III'  was  that  the  N2 
concentration  was  not  constant,  but  exhibited  a  marked  slope. 
There  are  two  possible  explanations  for  this  phenomenon. 
The  first  is  based  on  "regional  inhomogeneity , "  that  is, 
the  concentration  of  N~  in  the  alveoli  near  the  top  of  the 
lung  is  less  than  at  the  bottom  because  the  upper  regions 
are  initially  more  expanded  than  the  lower  due  to  the 
difference  in  pleural  pressure.   Thus,  less  oxygen  may 
enter  the  lower  regions  during  inspiration.   This  also 
requires  that  the  air  from  the  upper  regions  is  expired 
first,  often  called  "sequential  emptying."   There  is  evidence 
supporting  regional  inhomogeneity;  however,  sequential 
emptying  is  quite  controversial.   The  second  explanation  is 
that  of  "stratified  inhomogeneity,"  that  is,  there  is  a 
concentration  gradient  between  the  inspired  0   and  the  N 
rich  gas  already  present,  and  there  is  no  a  priori  reason  to 
assume  this  gradient  will  not  exist  in  the  alveoli. 

The  simple  model  which  follows  was  developed  to  study 
the  gas  mixing  problem,  and  the  N~  washout  case  in  particular, 
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B.   THE  SIMPLE  MODEL 

The  basic  development  of  this  model  was  suggested 
by  Pedley  [Ref.  17] . 

It  has  been  noted  that  the  lung  is  a  very  complicated 
network  of  branched  tubes,  in  which  the  exact  velocity 
distribution  remains  unknown.   Reynolds  numbers  may  vary 
from  several  thousand  to  less  than  one.   Thus,  solution  of 
the  full  diffusion  equation,  for  which  the  fluid  velocity 
field  is  required,  is  impossible,  and  considerable  simpli- 
fication is  necessary. 

The  lung  was  assumed  to  be  a  symmetric,  dichotomously 
branching  system,  in  which  all  pathways  have  the  same 
length.   The  lung  was  replaced  with  a  one-dimensional 
system  whose  cross-sectional  area  a(x)  increased  rapidly 
from  the  trachea  to  alveoli  as  shown  in  Figure  3  and 
consistent  with  Weibel's  [Ref.  21]  data.   The  area  was 
considered  constant  from  the  mouth  to  the  first  branching 

of  the  airways,  where  it  increased  greatly.   The  cross- 

2 
sectional  area  of  the  trachea,  an,  is  approximately  2.5  cm 

The  area  of  the  final  generation  of  airways  is  of  the  order 

4    2 
of  10   cm  .   The  position  where  the  area  began  to  increase 

was  taken  to  be  x  =  0.   The  distance  from  the  mouth  to 

x  =  0  is  about  38  cm,  while  the  remaining  distance  to  the 

final  generation  is  about  6.6  cm.   a(x)  was  assumed  to  be 

a  continuous  smooth  function  independent  of  time.   All 

volume  change  was  considered  to  take  place  in  the  alveoli. 
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The  N?  concentration  was  considered  to  be  uniform  across 
a  cross  section.   Thus,  the  average  concentration  c(x,t) 
was  a  function  only  of  distance  x  down  the  lung  and  of  time  t. 
The  velocity  u(x,t)  was  also  assumed  to  be  uniform  across 
a  cross  section.   These  assumptions  were  not  as  drastic  as 
they  first  appear.   The  trumpet-shaped'  tube  in  Figure  3  is 
actually  a  large  number  of  small  tubes  nearly  in  parallel, 
especially  at  the  alveolar  end.   The  main  effect  of  radial 
velocities,  velocity  profiles,  and  turbulence  is  to 
influence  the  relationship  of  concentration  to  the  longi- 
tudinal velocity.   Compensation  for  this  omission  was 
introduced  in  the  diffusivity. 

Previous  theories  by  Cuming,  Crank,  Horsfield,  and 
Parker  [Ref.  22]  and  La  Force  and  Lewis  [Ref.  23]  have  used 
similar  models,  but  greatly  oversimplified  the  actual  gas 
mixing  process.   They  ignored  the  fact  that  mixing  takes 
place  everywhere  simultaneously  with  convection,  especially 
in  the  upper  airways.   They  assumed  a  period  of  convection, 
in  which  the  interface  between  pure  0?   and  mixed  air  remained 
flat,  followed  by  a  period  of  diffusion,  followed  by  another 
period  of  convection.   Secondly,  they  assumed  that  the 
diffusion  was  governed  at  all  times  by  the  molecular  diffu- 
sivity D   , .   In  reality,  in  both  laminar  and  turbulent  flow, 
J      mol  J  '  * 

the  existence  of  a  velocity  profile  or  turbulent  eddy  ensures 
that  longitudinal  mixing  is  much  greater  than  can  be  described 
by  pure  molecular  diffusivity.   Taylor  [Refs.  24  and  25]  found 
that  after  a  certain  time  the  average  concentration  could  be 


21 


calculated  from  a  knowledge  of  the  average  longitudinal 

velocity,  on  the  basis  of  an  effective  longitudinal 

diffusivity  much  greater  than  D   , . 
7       fa  mol 

Flow  in  the  trachea  was  assumed  to  be  turbulent  and 
an  effective  longitudinal  diffusivity  (D~)  was  calculated 

using  Taylor's  relation.   For  a  tracheal  Reynolds  number  of 

2 
5000,  Dn  was  approximately  260  cm  /sec.   In  the  terminal 

regions  of  the  lung,  the  mixing  is  governed  almost  entirely 

2 
by  the  molecular  diffusivity  with  a  value  of  about  0.2  cm  / 

sec.   In  between  these  values,  the  flow  undergoes  large 

changes,  about  which  little  is  known  specifically.   D(x) 

was  considered  to  vary  smoothly  between  Dn  and  D   ,  from 

1  '  0      mol 

x  =  0  to  the  alveolar  end.   Therefore,  the  important 
feature  of  greatest  mixing  at  the  highest  velocities  was 
introduced. 

C.   ANALYSIS 

The  governing  equation  for  diffusion  in  the  one- 
dimensional  model  was: 


8c      3c    3_  rn  3c_.  +  D  9a  3c 
3t    U  3x  "  3x  L   3xJ    a  3x  3x 


where  the  last  term  accounts  for  the  change  in  cross 
sectional  area. 

From  continuity: 


U0a0 
u(x,t)  =  -f-i 
*  J        a(xj 
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where  un  was  the  velocity  in  the  trachea,  and  a^  was  the 
cross-sectional  area  of  the  trachea.   Inspiratory  and 
expiratory  flow  rates  were  considered  uniform  and  to  be  equal 
and  opposite. 

The  boundary  conditions  on  c(x,t)  were  as  follows. 
Both  c  and  dc/dx   were  continuous  at  x  =  0.   At  the  alveolar 
end,  for  N-  there  was  no  flux  across  the  end,  dc/dx   =  0. 
At  the  mouth,  c  =  0  for  inspiration  (pure  0~)  and  dc/dx   =    0 
for  expiration  (diffusion  negligible  compared  with 
convection) . 

The  initial  condition  for  inspiration  was  uniform 
concentration  everywhere.   The  initial  condition  for  expira- 
tion was  the  distribution  of  concentration  at  the  end  of 
inspiration. 

In  the  governing  equation  for  diffusion,  the  concentration 
was  non-dimensionalized  with  respect  to  its  initial  value, 
and  all  other  quantities  with  respect  to  their  values  in 

n  tun 

the  trachea  (C  =  —  ,D  =  ^-,X  =  £,A  =  —  ,T  =  -^-)  . 

c0       U0  *       a0 

The  following  equation  was  obtained: 

|C  AN  ♦  [±N  .  |_  (Afi)]  |C  .  AB  ifc 

where  the  minus  sign  denotes  inspiration,  the  plus  sign  for 
expiration,  and  N  =  un&/Dn,  the  effective  Peclet  number  in 
the  trachea. 
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The  new  boundary  and  initial  conditions  were 

C(X,0)  =  1 


for  inspiration 


C(-5,T)  =  0 

|£  (1-6, T)  =  0  ' 
for  expiration: 

|§  (-5.T)  =  |£  (1-6.T)  =  0 

where  6  =  (a  ,  /an)   ,  and  ensured  that  for  reasonable 
v  alv   0    ' 

Peclet  numbers,  the  values  of  A  and  D  were  close  to  their 
actual  values  (6  Z    0.02),  at  X  =  1-6. 

The  particular  forms  chosen  for  A(X)  and  D(X)  were: 

-5  <  X  <  0  A  =  D  =  1 

0  <  X  <  (1-6)        A  =  (1-X)"2 

1-D 

0  <  X  <  (1-6)         D  =  D   ,  +     , 
—   —  v    J  mol      .2 

_  2 
The  function  of  A  =  (1-X)   '  was  consistent  with  Weibel's 

data. 

D.   SOLUTION 

Solution  of  the  non-dimensional  equation  was  performed 
numerically  with  the  use  of  finite  difference  approximations: 
9C/3T  by  forward  difference,  and  8C/8X,  d2C/dX2   by  central 
difference.   The  resulting  equation  was  solved  on  the  digital 
computer  using  sufficiently  small  time  and  distance  steps  to 
maintain  stability  and  accuracy.   The  concentration  at  any 
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time  T  was  calculated  from  data  for  the  previous  time  step. 
Solutions  were  generated  for  four  values  of  the  effective 
Peclet  number  in  the  trachea.   Inspiration  time  before 
expiration  was  also  varied. 

E.   RESULTS  AND  DISCUSSION 

The  simplified  model  of  the  lung  produced  results 
similar  to  those  described  by  Farhi  [Ref.  16]  and  Pedley 
[Ref.  17],  with  a  minor  variation.   In  addition,  several 
interesting  results  were  observed. 

The  effective  Peclet  number  in  the  trachea  was  found 
to  have  an  effect  that  varied  with  distance  down  the  model, 
as  shown  in  Figure  4  and  by  data  in  Appendix  A.   In  the 
trachea,  an  increase  in  the  Peclet  number  caused  a  decrease 
in  the  N~  concentration  at  a  particular  position.   This 
effect  was  reversed  for  all  X  greater  than  zero,  in  the 
actual  branching.   That  is,  an  increase  in  Peclet  number 
increased  the  N~  concentration  at  a  particular  position. 
In  addition,  this  effect  tended  to  be  felt  just  upstream  of 
X  =  0.   The  general  result  was  that  as  the  Peclet  number 
decreased,  the  further  down  the  lung  the  diffusion  progressed, 
This  indicates  that  a  reduced  Peclet  number  corresponded  to 
relative  increase  in  the  diffusivity  with  respect  to  the 
velocity. 

The  effect  of  the  duration  of  inspiration  was  also 
dependent  on  distance  down  the  model.   In  the  trachea,  the 
concentration  rapidly  decreased  with  time,  while  in  the 
lung  the  process  was  much  slower,  as  shown  in  Figure  5. 
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N   concentration  in  the  model  during  expiration  also 
varied  considerably  from  the  trachea  to  the  lung.   As 
expiration  time  increased,  the  N?  concentration  slowly 
increased  linearly  along  the  trachea,  then  rapidly  increased 
with  distance  to  the  alveoli  (Figures  6,  7,  8).   Increasing 
expiration  time  also  produced  an  increase  in  the  N?  concen- 
tration at  a  particular  position,  with  the  exception  of  the 
last  few  airways  before  the  alveoli. 

The  N~  washout  curves  obtained  were  similar  to  those 

expected  with  the  exception  that  the  sharper  changes  in 

slope  were  not  detectable  between  regions  II  and  III,  as 

shown  in  Figures  9,  10,  and  11.   In  fact,  for  the  longest 

inspiration  (T  =  96) ,  the  curve  for  Peclet  number  of  12 

was  essentially  linear.   For  this  reason  it  was  very 

difficult  to  examine  the  relationships  between  V~n ,  V~, 

and  V-   as  discussed  in  the  introduction.   V_.  was  55 
in  D 

(non-dimensional),  and  for  Figure  9,  V~n  =  12,  and  the 
relationship  appeared  to  hold.   For  Figure  10,  V.   =  V~n  =  48, 
but  it  was  difficult  to  determine  exactly  where  the  alveolar 
plateau  began.   For  Figure  11,  Vn  =  96  and  the  alveolar 
plateau  should  have  begun  at  Vn  =  55.   It  obviously  did  not 
for  either  case. 

The  effect  of  Peclet  number  was  to  decrease  the 
concentration  at  a  particular  time  with  an  increase  in 
Peclet  number  for  the  washout  curves. 
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F.   CONCLUSIONS 

The  results  of  the  numerical  solution  for  the 
theoretical  simple  model  of  gas  mixing  in  the  lung  justified 
the  following  conclusions: 

1.  The  model  effectively  simulated  gas  mixing  in  the 
lung  with  simultaneous  diffusion  and  convection  effects. 

2.  The  model  is  very  flexible  and  many  parameters  may 

be  easily  varied  to  more  closely  approximate  the  experimental 
data  for  the  washout  curves. 

3.  The  flexibility  of  the  model  extends  to  many  other 
similar  experiments  in  the  field.   The  ease  with  which 
initial  conditions,  boundary  conditions,  and  the  various 
parameters  may  be  varied  and  the  rapidity  of  computer 
solution  provide  a  fast,  accurate  tool  for  future  research 
in  the  field. 
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III.   EXPERIMENTAL  STUDY  OF  FLOW 
IN  A  THIN-WALLED  COLLAPSIBLE  TUBE 

A.   INTRODUCTION 

It  was  noted  earlier,  in  the  survey  section  of  this 
paper,  that  flow  in  the  small  capillaries  of  the  lung  is 
analogous  to  flow  in  a  thin-walled  collapsible  tube.   Flow 
in  this  type  of  tube  is  dependent  on  the  surrounding  pressure 
as  well  as  the  upstream  and  downstream  pressures. 

The  force  of  gravity  exerted  on  the  lung  tends  to 
divide  it  into  three  basic  zones.   In  zone  I  at  the  top  of 
the  lung,  the  alveolar  pressure  PA  is  greater  than  the  up- 
stream static  pressure  P, .   In  zone  II,  near  the  middle  of 
the  lung,  P,  is  greater  than  P.  which  is  in  turn  greater 
than  the  downstream  static  pressure  P?.   In  the  lowest 
part  of  the  lung,  zone  III,  P.  is  less  than  both  P,  and 
P~.   The  pressure- flow  relationships  in  these  three  zones 
have  been  experimentally  studied  by  Conrad  [Ref.  8]  and 
Holt  [Ref.  9]  and  many  others.   The  results  have  all 
indicated  that  the  flow  is  pressure  dependent  as  follows. 
In  zone  I  the  tube  was  collapsed  and  there  was  no  flow.   In 
zone  II  the  flow  was  proportional  to  P,  -  P.,  and  the  tube 
was  observed  to  neck  near  the  downstream  end.   In  zone  III 
the  flow  was  unobstructed  and  proportional  to  P,  -  P~. 
These  results  were  obtained  using  a  latex  rubber  tube  as  a 
model  surrounded  by  an  airtight  box.   The  fluid  involved 
was  water.   The  results  were  based  on  a  single  length-to- 
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diameter  (L/D)  ratio  and  a  single  tube  tension  in  each 
study. 

This  experiment  was  conducted  using  air  as  the  fluid 

in  order  to  verify  the  previous  results  and  to  study  the 

effects  of  variations  in  L/D  ratio  and  tension  in  the  tube 
wall. 

B.   EXPERIMENTAL  APPARATUS 

The  experimental  apparatus  used  in  this  investigation 
was  designed  and  constructed  similar  to  the  previous  models 
of  Conrad  [Ref.  8]  and  Holt  [Ref.  9].   Alterations  were 
introduced  where  necessary  to  simplify  the  apparatus  and 
adapt  it  to  air  flow.   The  system  is  shown  in  Figures  12 
and  13.   The  details  of  the  equipment  were  as  follows. 

An  airtight  box  was  constructed  with  dimensions  14"  by 
10"  by  10"  using  1/2"  Plexiglas  sheet.   The  inlet  to  the 
box  was  a  decreasing  radius  nozzle  constructed  from  powdered 
aluminum- filled  epoxy  resin.   It  converged  smoothly  from  an 
outside  diameter  of  5-1/2"  to  an  inner  diameter  of  1" 
aluminum  pipe.   The  outside  diameter  of  the  pipe  was  1-1/4". 
The  exit  from  the  box  was  via  the  same  type  pipe.   The  mouths 
of  the  two  pipes  were  separated  by  10",  and  radiused  to 
ensure  smooth  flow.   Variations  in  the  L/D  ratio  were 
attained  using  several  lengths  of  Plexiglas  tubing  adjusted 
over  the  aluminum  pipe.   The  inside  diameter  of  the 
Plexiglas  tubing  was  1-1/4",  and  the  outside  diameter  was 
1-1/2".   The  mouth  of  each  tube  was  also  radiused.   A  1" 


29 


gate  valve  was  mounted  between  the  exit  and  the  vacuum 
source  in  order  to  adjust  the  pressure.   Pressure  in  the  box 
was  regulated  by  a  3/4"  bleed  from  the  exit  pipe  (taken 
upstream  of  the  1"  gate  valve)  equipped  with  a  3/4"  gate 
valve.   Static  pressure  ports  were  located  at  the  inlet, 
exit,  and  top  of  the  box.   A  pitot  tube  was  inserted  in 
the  exit  flow.   The  pressure  ports  were  connected  to  an 
inclined  water  manometer.   The  suction  source  was  a  Singer 
vacuum  cleaner,  Model  D-5.   The  collapsible  tube  was  a  thin 
rubber  prophylactic  with  a  limp  length  and  diameter  of  7" 
and  1-3/8",  respectively.   It  was  stretched  over  the 
Plexiglas  tube  at  each  end,  and  fixed  by  rubber  O-rings  and 
bands . 

C.   EXPERIMENTAL  PROCEDURE 

The  basic  experimental  procedure  consisted  of  variation 
of  tube  tension  for  each  of  three  L/D  ratios. 

The  1"  gate  valve  was  opened  to  adjust  the  inlet  pressure 

P..  to  a  value  of  2.0  in  H?0  below  atmospheric  pressure,  and 

maintained  at  this  level  for  all  investigations.   One  of 

three  L/D  ratios  was  then  chosen  and  the  tube  tension  was 

adjusted  for  each  of  three  values.   A  new  L/D  ratio  was 

chosen  and  the  procedure  was  repeated.   The  pressure  in  the 

box  P.  was  varied  from  5.0  in  H?0  below  atmospheric  up  to 

collapse  of  the  tube.   As  a  tube  was  stretched,  its  diameter 

decreased,  and  the  tension  of  the  tube  was  quantified  by 

the  pressure  difference  P,  -  P   necessary  to  maintain 

1    A 
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parallel  flow  in  the  tube  (tube  walls  parallel  to  the  pipe 
walls  at  the  same  diameter) . 

D.   RESULTS  AND  DISCUSSION 

The  experimental  study  of  air  flow  in  a  collapsible  tube 
produced  significant  results.   The  data  collected  are  tabu- 
lated in  Appendix  B  and  graphically  shown  in  Figures  14,  15, 
and  16.   In  general,  the  results  obtained  were  similar  to 
those  of  previous  work  with  water,  and  additional  new  infor- 
mation was  observed. 

In  zone  I  the  flow  was  discovered  to  differ  from  zero, 
depending  on  the  initial  tension  of  the  tube.   In  the  low 
tension  case,  flow  was  found  to  exist  for  P.  greater  than 
P1  before  collapse  of  the  tube.   Generally,  the  effect  of 
decreasing  tension  was  to  increase  the  pressure  P.  attained 
before  collapse.   In  zone  II  the  slope  of  the  dynamic  pres- 
sure curve  (directly  related  to  the  mass  flow)  rapidly  in- 
creased as  P.  passed  P?,  and  leveled  as  P.  approached  (and 
passed  in  some  cases)  P, .   In  zone  III  the  major  deviation 
from  previous  results  occurred.   As  the  difference  P,  -  P- 
was  maintained  as  a  constant  in  each  case,  the  dynamic  pres- 
sure should  have  remained  constant  for  all  P.  less  than  P,. 
The  results  of  this  investigation  indicated  that  the  flow 
was  constant  until  parallel  flow  was  attained.   However, 
further  decrease  in  P.  caused  the  cross  section  to  balloon, 
accompanied  by  a  drop  in  the  dynamic  pressure. 
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As  depicted  in  Figures  14,  15,  and  16,  the  effect  of 
varying  L/D  ratios  was  to  increase  the  P.  possible  before 
collapse,  with  decreasing  L/D. 

Three  qualitative  phenomena  were  observed  during 
experimentation.   First,  just  prior  to  collapse,  the  tube 
walls  developed  longitudinal  striations  near  the  center  of 
the  tube.   Second,  collapse  was  a  highly  non-steady  phenome- 
non.  It  was  not  complete,  steady  collapse,  but  a  rapid 
oscillatory  motion  of  the  tube  walls  from  complete  collapse 
to  the  striated  state  described  above.   Finally,  the  col- 
lapse was  observed  to  take  place  at  a  position  approximately 
1/5  of  the  tube  length  from  the  downstream  end  for  each 
value  of  L/D.   The  "ballooning"  effect  was  also  observed  to 
begin  at  the  downstream  end. 

E.   CONCLUSIONS 

The  results  of  the  investigation  of  air  flow  in  a  thin- 
walled  collapsible  tube  justified  the  following  conclusions: 

1.  The  experimental  apparatus  used  was  an  accurate 
model  of  flow  in  a  collapsible  tube  utilizing  new  parameter 
variations . 

2.  The  length-to-diameter  ratio  of  the  flexible  tube 
was  found  to  have  a  significant  effect  on  pressure- flow 
relationships . 

3.  The  longitudinal  tension  of  the  flexible  tube  was 
also  found  to  have  a  very  significant  effect  on  the 
pressure-flow  relationships. 
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FIGURE  17.   "PARALLEL"  FLOW  IN  COLLAPSIBLE  TUBE 
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FIGURE  18.   "BALLOON"  FLOW  IN  COLLAPSIBLE  TUBE 


FIGURE  19.   STRIATIONS  PRIOR  TO  COLLAPSE 
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FIGURE  20.   UNSTEADY  COLLAPSE 
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APPENDIX  A 
NUMERICAL  SOLUTION  DATA 


N~  Concentration  during  Inspiration,  T.   =4.8 

2  r       »   in 


Distance 

(x/£) 

N 

=    6 

-5.0 

.000 

-4.0 

.000 

-3.0 

.006 

-2.0 

.050 

-1.4 

.130 

-1.0 

.  216 

-0.4 

.389 

0.0 

.545 

0.2 

1 

.000 

0.6 

1 

.000 

0.98 

1 

.000 

in 

Concentration  (c/cn) 

N  =  8  U    N  =  10  N  =  12 

.000  .000  .000 

.000  .000  .000 

.002  .000  .000 

.030  .017  .008 

.103  .087  .062 

.193  .172  .154 

.387  .387  .386 

.549  .554  .559 

1.000  1.000  1.000 

1.000  1.000  1.000 

1.000  1.000  1.000 
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N?  Concentration  during  Inspiration,  N  =  10 


Distance  Concentration  (c/cn) 


(x/£) 

T.   =  4.8 
in 

T.   =  12 
in 

T.   =  48 
in 

T.   =  96 
in 

5.0 

.000 

.000 

.000 

.000 

4.0 

.000 

.000 

.000 

.000 

3.0 

.000 

.000 

.000 

.000 

2.0 

0.17 

.000 

.000 

.000 

1.4 

.081 

.000 

.000 

.000 

1.0 

.172 

.000 

.000 

.000 

0.4 

.387 

.000 

.000 

.000 

0.0 

.554 

.000 

.000 

.000 

0.1 

.999 

.870 

.301 

.129 

0.2 

1.000 

.998 

.641 

.321 

0.4 

1.000 

1.000 

.970 

.746 

0.6 

1.000 

1.000 

.999 

.966 

0.8 

1.000 

1.000 

1.000 

.999 

0.98 

1.000 

1.000 

1.000 

1.000 
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N?  Concentration  during  Expiration,  N  =  10 


Distance  Concentration  (c/cn) 


(x/JO 

T.  =12 
in 

T.  = 
in 

=  48 

T.  =96 
in 

T   =12 
ex 

T   =12 
ex 

T   =48 
ex 

T   =12 
ex 

T   =48 
ex 

T   =96 
ex 

-5.0 

.532 

.170 

.586 

.073 

.344 

.642 

-4.0 

.559 

.185 

.593 

.080 

.351 

.646 

-3.0 

.586 

.201 

.601 

.087 

.358 

.646 

-2.0 

.610 

.216 

.609 

.095 

.365 

.656 

-1,0 

.632 

.231 

.616 

.102 

.373 

.661 

0.0 

.652 

.246 

.623 

.110 

.380 

.666 

0.1 

.853 

.439 

.720 

.216 

.479 

.728 

0.2 

.984 

.719 

.859 

.418 

.650 

.828 

0.3 

.999 

.898 

.945 

.628 

.799 

.904 

0.4 

1.000 

.973 

.982 

.800 

.894 

.949 

0.6 

1.000 

.999 

.998 

.969 

.979 

.988 

0.8 

1.000 

.1.000 

1.000 

.999 

.998 

.998 

0.98 

1.000 

1.000 

1.000 

1.000 

1.000 

.999 
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Washout  Data,  T.   =  12 
'   in 


T 
ex 

uQt/£) 

2. 

4 

4. 

8 

7. 

,2 

9. 

,6 

12. 

,0 

14. 

,4 

16. 

,8 

19. 

,2 

21, 

,6 

24, 

,0 

Concentration  (c/cn) 
N  =  6  N  =  12 


000  .000 

101  .061 

317  .298 

453  .438 

540  .530 

603  .597 

651  .650 

689  .692 

721  .726 

748  .755 
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Washout  Data,  T.   =  48 
*   in 


ex 
(uQt/£) 

2.4 

4.8 

6.0 
12.0 
18.0 
24.0 
30.0 
36.0 
42.0 
48.0 
54.0 
60.0 
66.0 
72.0 


Concentration  (c/cn) 


N  =  6 

.000 
.030 
.067 
.207 
.304 
.384 
.453 
.513 
.566 
.612 
.653 
.688 
.720 
.748 


N  =  12 

.000 
.013 
.042 
.157 
.247 
.328 
.400 
.466 
.525 
.578 
.626 
.668 
.705 
.738 
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Washout  Data,  T.   =96 

in 


T 
ex 

(u0t/£) 

2. 

4 

4. 

8 

6. 

0 

9. 

8 

12. 

0 

24. 

0 

36. 

0 

48. 

0 

60. 

,0 

72. 

,0 

84. 

,0 

96. 

,0 

108, 

,0 

120. 

.0 

Concentration  (c/cn) 
N  =  6  .  n  =  12 


000  .000 

010  .002 

034  .016 

086  .046 

110  .061 

224  .144 

325  .231 

417  .320 

499  .407 

571  .488 

634  .562 

688  .629 

735  .687 

774  .737 
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APPENDIX  B 
THIN-WALLED  COLLAPSIBLE  TUBE  DATA 


All  pressures  are  inches  of  water  below  atmospheric 
pressure.   P1  is  2.0  for  all  cases. 


=  2.3 


L/D   = 

5.5,  AP 
P 

P2 

PA 

2.85 

5.0 

2.85 

4.5 

2.85 

4.3 

2.85 

4.0 

2.85 

3.5 

2.85 

3.0 

2.85 

2.8 

2.85 

2.6 

L/D   = 

5.5,  AP 
P 

2.80 

5.0 

2.80 

- 

4.5 

2.80 

4.0 

2.80 

3.9 

2.80 

3.5 

2.80 

3.0 

2.80 

2.6 

2.80 

2.4 

2.80 

2.25 

2.80 

2.2 

=  1.9 


65 

60 
60 
55 
50 
45 
45 
40 


65 
60 
55 
55 
55 
50 
45 
40 
40 
35 
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L/D  =  5.5,  AP   =1.5 
P 


r2 

A 

2.75 

5.0 

2.75 

4.5 

2.75 

4.0 

2.75 

3.5 

2.75 

3.0 

2.75 

2.75 

2.75 

2.2 

2.75 

2.0 

2.75 

1.9 

q 

65 
65 
60 
55 
55 
50 
45 
40 
40 


L/D  =  4.0,  AP   =  2.3 
P 


2.80 

5.0 

2.80 

4.5 

2.80 

4.3 

2.80 

4.0 

2.80 

3.5 

2.80 

3.0 

60 
55 
50 
45 
45 
45 


55 


L/D  =  4.0,  AP   =  1.9 
P 


r2 

A 

2.80 

5.0 

2.80 

4.5 

2.80 

4.0 

2.80 

3.9 

2.80 

3.5 

2.80 

3.0 

2.80 

2.8 

2.80 

2.5 

2.80 

2.1 

q 

60 

50 

45 

45 

45 

.40 

35 

,35 

,35 


L/D  =  4.0,  AP   =  1.5 


2.75 

5.0 

2.75 

4.5 

2.75 

4.3 

2.75 

4.0 

2.75 

3.5 

2.75 

3.0 

2.75 

2.75 

2.75 

2.5 

2.75 

2.2 

2.75 

2.0 

2.75 

2.0 

2.75 

2.1 

50 
45 
40 
40 
40 
40 
40 
40 
35 
30 
30 
30 
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L/D  =  2.5,  AP   =  2.3 
P 


r2 

A 

2.75 

5.0 

2.75 

4.5 

2.75 

4.3 

2.75 

4.0 

2.75 

3.5 

2.75 

3.0 

2.75 

2.75 

2.75 

2.2 

2.75 

2.0 

q 

50 
45 
40 
40 
40 
40 
40 
35 
35 


L/D  =  2.5,  AP   =1.9 


2.70 

5.0 

2.70 

4.5 

2.70 

4.0 

2.70 

3.9 

2.70 

3.5 

2.70 

3.0 

2.70 

2.7 

2.70 

2.4 

2.70 

2.0 

2.70 

1.8 

2.70 

1.6 

50 
45 
40 
40 
40 
40 
40 
40 
35 
35 
35 
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L/D   =    2.5,    AP      =1.5 
P 


2 

XA 

2.65 

5.0 

2.65 

4.5 

2.65 

4.0 

2.65 

3.5 

2.65 

3.0 

2.65 

2.6 

2.65 

2.0 

2.65 

1.8 

2.65 

1.5 

2.65 

1.6 

q 

.50 
.45 
.40 
.40 
.35 
.35 
.35 
.30 
.30 
.30 
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COMPUTER  NUMERICAL  SOLUTICN  FCR  GAS  MIXING 
IN  A  SIMPLE  MCDEL  OF  THE  LUNG 


DIMENSION  CNW(124)  ,CTN(124) 

NP=10 

WRITE  (6,250)  NP 
25C  FCR^AT  (5X,'PECLET  NO.  IS*  ,3X , 13  ,// ) 

T  =  C.O 

CELX=0.2 

JE=124 

NM=2 

JEi=JE-l 

DC  10  1=1,4000 

X=-5.0 

D=1.0 

F=-1.0 

A=1.0 

DELX=.2 

DELT=( .3*(DELX**2) ) /D 

IF  (I  .GT.l)  GO  TO  111 

WRITE  (6,150) 
150  FORMAT  (5X,'T  EQUALS    0.0',/) 

WRITE  (6,300) 

SET  THE  INITIAL  CONDITIONS 

DC  110  K=1,JE 
CTN(K)=1.0 

IF  (X.LT.O.)  GO  TO  30 

DELX=.01 
3C  X=X+DELX 

WPITE  (6,310)  X,CTN(K) 
11C  CONTINUE 

X=-5.0 

DELX=.2 
111  CTN(l)=C.O 

T=T+DELT 

MCC=( 1/500  )*500 

IF  (MOD.EQ.I)  GO  TO  50 

GO  TO  70 
50  WRITE  (6,200)  T 
20C  FORMAT  (//,5X,'T  E QUALS' ,3X ,F8 . 4  ,  / ) 

WRITE  (6,300) 
30C  FCRMAT  (5X,'X  EQUALS ', 5X, • C  EQUALS',/) 
70  DC  100  J=2,JE1 

CALCULATE  CONCENTRATION  FROM  MOUTH  TO  ALVECLI 
FOR  EACH  TIME  STEP 

11  IF  (X.LT.O.)  GO  TO  40 

DELX=.0  1 

DELT=(.3*(CELX**2) )/D 

A=( l-X)**(-2) 

D=.2+.8*( U-X)**(2*NM) ) 

FB=( .4+1.6*(l-X)**(2*NM) ) / ( NP* ( 1-X ) )-l/ A 

F= (NM/N P )*<-!. 6) *(  1-X) **(2*NM-1)+FB 

/a.  r*     y  —  y  j.nt  i  v 

G=(D/NF)*(CTN(J+1 )-2.*CTN(J)+CTN(J-l) )/ (DELX**2) 

H=(CTN( J  +  l )-CTN( J-l)  )/(2.*DELX) 

CNW(J)=CTN(J)+DELT*(G+(F*H) ) 

IF  (MOD.EQ.I )  GO  TC  55 

GO  TO  100 
55  WRITE  (6,310)  X,CNW(J) 
31C  FORMAT  (  5X ,F 8.  4, 3X ,F 12 . 7) 
100  CCNTINUE 

X=X+DELX 

IMPOSE     B.C.       DC/DX=0       AT    ALVOELAR    END 

CNW(JE )=4./22.*( 9*CNW( JE 1 ) -  18 . / 4 . *CNW ( JE-2  MCNW  ( J  E-3  ) ) 
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DC    7    L=2tJE 

C7ML)  =CNW(L) 
7    CONTINUE 

IF     (MOD. EG. I)    GO    TO    60 

GC  TO  10 
60  WRITE  ( 6,310)  X,CTN( JE  ) 
10  CONTINUE 

STCP  INSPIRATION,  START  EXPIRATION 

CTN(l)=C..O 

WRITE     (6,700) 
700    FORMAT     (//,5X,'STCP     INSPIRATION,     START    EXPIRATION',//) 

DO    15    1=1,6000 

DELX=.2 

X=-5.0 
D-1.0 

A=1.0 

F  =  1.0 

DELT=( ,3*(DELX**2)  ) /D 

T=T+DELT 

IF    (I .GT.500)    GO    TO    17 

MAD=( I/20)*20 

IF    (MAD.NE.I )     GO    TO    17 

WRITE    (6,200)     T 

WRITE     (6,300) 

WRITE     (6,310)     XfCTNCU 
17    MOD=( I/500)*50O 

IF     (MOD.EQ.I )    GO    TO    205 

GC  TO  75 
205  WRITE  (6,200)  T 

WRITE  (6,300) 

WRITE  (6,310)  XtCTMl) 
75  DO  105  J=2,JE1 

CALCULATE  CONC ENTRATI CN  FROM  MOUTH  TO  ALVECLI 
FOR  EACH  TIME  STEP 

IF  (X.LT.O.)  GO  TO  45 

DELX=.01 

DELT=( .3*(DELX**2)  )/D 

A  =  (  l-X)*=M-2) 

D=.2+.6*< (  l-X)**(2*NM)  ) 

FE=( .4+l.6*(l-X)**(2*NM) )/ (NP*(  1-X ) )+l/A 

F=(NM/NP)*(-1.6) *(  1-X)**(2*NM-1 )+FB 
45  X=X+DELX 

G=(D/NP)*(CTN(J+1)-2.*CTN(J )+CTN(J-l) )/ (DELX**2) 

H=(CTM  J+1)-CTN(J-1)  )/(2.*DELX) 

CNW(J)=CTN(J)+DELT*(G+(F*H)  ) 

IF  (MOD. EC. I)  GO  TO  90 

GC  TO  1C5 
90  WRITE  (6,310)  X,CNW(J) 
1C5  CCNTINUE 

X=X+DELX 

IMPOSE  B.C.   DC/DX=0   AT  ALVOELAR  END 

CNW(JE)  =  4./22.=M9*CNW(  JED-  18 . /4 .*CNW( JE-2 ) +CNW( JE-3) ) 

IMPOSE  B.C.   DC/DX=0   AT  MOUTH 

CTN(1 )=4./22.*(9 .*CNW(2)-18 ./4. *CNW ( 3) +CNW ( A )  ) 

DO  9  L=2,JE 

CTN(L)=CNW(L) 

9  CCNTINUE 

IF  (MOD.EQ.I )  GO  TC  65 

GC  TO  15 
65  WRITE  (6,310)  X,CTN(JE) 
15  CCNTINUE 

STCP 

END 
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